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Abstract 

We outline how relic neutrino asymmetries may be generated in the early universe 
via active-sterile neutrino oscillations. We discuss possible consequences for big bang 
nucleosynthesis, within the context of a particular 4-neutrino model. 



1 Introduction 



The implications of neutrino masses in cosmology and astrophysics are numer- 
ous. One fascinating effect, is that of relic neutrino asymmetries, which may 
be produced if there is mixing between active and sterile neutrino species [1]. 

Neutrino asymmetries are interesting since they affect big bang nucleosynthe- 
sis (BBN), altering the primordial helium yield. BBN is particularly sensitive 
to an asymmetry between the v e and V e . Additionally, since the matter poten- 
tial arising from forward scattering is proportional to the asymmetry, a large 
asymmetry will suppress active-sterile oscillation modes which could otherwise 
equilibrate sterile neutrino species around the time of BBN. 



* Invited talk at the NUFACT'01 workshop in Tsukuba, Japan 24-30 May 2001. 
Email address: n.bell@physics.unimelb.edu.au (Nicole F. Bell). 



Preprint submitted to Elsevier Science 



7 February 2008 



2 Relic neutrino asymmetries 



First consider two-flavour mixing between an active neutrino, u a , and a sterile 
neutrino, v s . For each neutrino flavour we define an asymmetry 

L a = (1) 

n 7 



where n v and n 7 are the neutrino and photon number densities, respectively. 
The neutrino momentum distributions, N(p), are such that 

n = J N{p)dp, (2) 



and in thermal equilibrium are just given by Fermi-Dirac distributions 

•V,(/') = ~ (3) 

2vr 2 1 + eX p (J + /i 



Here T is the temperature and p, = /x/T, with /i being the neutrino chemical 
potential. In equilibrium, L a is related to // as per 



L ~ 



1 



24C(3) 



7i 2 ((l u - ftp) - 6(Pl - In 2 + (p* - fil) . (4) 



We wish to track L a as the universe evolves, with the period of interest being 
the epoch between T ~ lOOMeV and T ~ lMeV. The evolution of the neutrino 
ensemble will be controlled by (i) oscillations, (ii) collisions of the neutrinos 
with other particles, resulting in decoherence of the oscillations, and (iii) the 
expansion of the universe, which redshifts energies and momenta. 

The matter affected mixing angle, 9 m , is 

• 2 0Q sin 2 26 

sin W m - g . n2 26o + {b±a _ cog 20o)2 , W 



where 6q is the mixing angle in vacuum. The term (b ± a) arises from the 
effective matter potential V = (b ± a)Sm 2 /2p, where -/+ corresponds to vjv 
and 5m 2 is the mass squared difference of the two neutrino mass eigenstates. 
The functions a and b are given by 



4V2((3)G F T 3 L^p 
7r 2 5m 2 
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4v/2C(3)GfTM qP 2 

b{p) = ' (6) 

where A e (A TiM ) ~ 17(4.9) [2]. The function L<°) is defined as 

= L a + L T + L M + L e + 77, (7) 



where 77 ~ 10~ 10 is due to the baryon asymmetry. Note that the mixing angle 
is dependent on the size of the asymmetry, which makes the evolution of the 
asymmetry non-linear. 

There are two distinct phases in the growth of a neutrino asymmetry. The 
first is a collision dominated stage, which occurs at high temperature. The 
total collision rate of a v e (z/ M or v T ) of momentum p with particles in the 
background plasma is approximately 

r(p)^4.0(2.9)^T 5 . (8) 



These collisions tend to destroy the coherent superposition of flavour states. 
In this regime, it can be shown that the growth of the neutrino asymmetry is 
governed by the equation 



dL a 71-2 f, ( AT + AT +\ s 2 Fa(c — b) 

~df ~ 2C(3)T 3 J P \ a ~ s ) [x + (c - b + a) 2 ][x + (c - b - a) 2 ]' ^ 

where c = cos 29 , s = sin 29 , 



x 



pr(p) 

Am 2 



1 2 



, N + = ^{N V + N 17 ). (10) 



Eq.(9) can be derived either from a heuristic picture [1] that views collisions 
as measurement-like interactions, or can be extracted from the exact Quan- 
tum Kinetic Equations [3]. It is valid when the collision rate is large and the 
evolution is adiabatic - see [3] for details. 

Note that dL a /dt oc a oc lS a \ which can lead to runaway positive feedback 
and thus rapid growth of L a . The mixing parameters for which this occurs are 
specified by \5m 2 \ ~ 10~ 4 eV 2 and sin 2 29 ~ 1CT 10 , and we require lighter of 
the two mass eigenstates to approximately coincide with the v s . Typically, the 
asymmetry will reach values L a ~ O(10~ 5 ) in this phase. Note that the sign 
of the asymmetry, which depends on the initial conditions, can be ambiguous 
[6]. 
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The second phase is that lower temperature. Here the collisions are less im- 
portant, and coherent MSW transitions drive the asymmetry to values L a ~ 
O(0. 1) [7] .Pi If, for example, the asymmetry is positive, the neutrino MSW 
resonance occurs at a much higher momentum value than that of the an- 
tineutrinos and may be neglected. The rate at which the asymmetry grows is 
then 



dT 



-X 



d /Pres 

dT \ T 



where X is the difference between the number of u a and u s at the resonance. 
Eq.(ll) quantifies the rate at which the resonance moves through the momen- 
tum distribution, and holds when transitions are adiabatic and the resonance 
width is small. 



3 Big bang nucleosynthesis and a four neutrino model 



We wish to determine the impact of such asymmetry production on the BBN 
light element abundances, and in particular, on the helium yield. In general, 
the results are somewhat model dependent (compare, say, [7] with [8].) 

The results are particularly sensitive to the size of the v e asymmetry, which 
directly affects the neutron/proton ratio through the reaction rates 



X(n — > p) ~ \{nu e — > pe ) + A(ne + — > pv e 
\(p — > n) ~ \(pe~ — > nv e ) + A(pz/ e — > ne + ). (12) 



The n/p ratio determines the helium mass fraction Yp via the equation 
dY 

_^ = _ A (n -> p)Y P + X(p - n)(2 - Y P ). (13) 
dt 



For typical values of 5m 2 , one must compute the growth of the asymmetry and 
the effect on BBN simultaneously, as the asymmetries will not have reached 
their final values before the BBN epoch commences. 

Momentum modes which are depleted due to the oscillations, are refilled by 

1 Ref.[4] disputed the production of large asymmetries. However, this issue has been 
resolved, with [4] shown to be in error due to the use of an approximation outside 
its region of validity [5] . 
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scattering and annihilation processes as per 

d ( -^actual \ T-i / \ / -^eq -^actual \ / , . \ 

Jt\—) =T{P) {1^-—)' (14) 



where N eq is the equilibrium distribution in Eq.(3), and N = A^ cq |- =0 . 

The example we consider is a particular "2+2" style 4-neutrino model, consist- 
ing of a maximally mixed — v T pair, separated by a mass gap from u e — v s . 
The oscillations modes most important in the evolution of the various flavour 
neutrino asymmetries will be the ones with the largest 5m 2 . Specifically, the 
— v s oscillation modes will create large and L T asymmetries, and the 
v^t — v e oscillation modes will transfer a small amount of this asymmetry to 
the electron neutrinos. See [8] for further details. 

Numerically integrating the appropriate generalisations of eq.(ll), we find that 
the magnitude of the final lepton asymmetries generated are 



\L„\/h= |L T |//i~0.16, 

\L e \/h~6.7 x KT 3 , (15) 

for 5m 2 ~ 4eV, where h = Tj^/T^. 

Integrating the BBN reaction rates, eq. (12), the direct effect of the v e asymme- 
try on the n/p ratio, and hence Yp, is determined. We obtain 5Yp ~ ±0.0023, 
corresponding to an effective change in the number of neutrinos of ~ 
±0.19. In addition, there is a small direct change to the energy density, leading 
overall to 



8NI — —0.3, positive asymmetry, 

5N° S ~ +0.1, negative asymmetry, (16) 
for 5m 2 ~ 3eV 2 . 

The value of N„ obtained in this way must then be compared with that 
calculated from the observed light abundances [9] . It is not yet clear whether 
the observed helium and deuterium abundances can be accommodated within 
the standard BBN scenario with = 3, and it is possible that a small v e 
asymmetry such as that discussed here would help to provide a better fit [10]. 
Recent measurement of the cosmic microwave background radiation (CMBR) 
anisotropy also provides a determination of the number of neutrinos in the 
early universe [11], although the constraints yet, not as strong as those 

arising from BBN.[3] However, future measurements are expected to pin down 

2 Note, however, that the effective number of neutrinos extracted from the CMBR 
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the effective neutrino number very precisely [12] and thus will be of great use 
in constraining models which feature the mixing of active and sterile neutrinos. 
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and BBN data are not necessarily the same parameter. The two determinations are 
complementary, since the BBN quantity depends on flavour through its sensitivity 
to the f e asymmetry, while the CMBR measures the total energy density residing 
in neutrinos (at a much later epoch of the universe.) 
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